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Abstract. During the last decade, the FORS 1 instrument of the ESO Very Large Telescope has been 
used to obtain low resolution circular polarized spectra for about 500 stars, with the aim of mea- 
suring their mean longitudinal magnetic fields. Magnetic field estimates were obtained by different 
authors, using different software tools. Several interesting detections were obtained at a 3 a level; 
some of them were eventually confirmed by follow-up investigations, some of them were not. This 
raises issues about the reliability of the stated uncertainties of some of the published field values. 
To investigate these problems, we have developed a semi-automatic procedure for magnetic field 
determination, which includes self-consistent checks for field detection reliability. We have applied 
our procedure to the full content of single star ("fast mode") circular spectropolarimetric measure- 
ments of the FORS 1 archive, and explored the details and interagreement of various methods for 
data reduction. We have finally produced a catalogue of FORS 1 longitudinal field measurements 
which includes about 1000 entries. Here we critically review the previously published FORS 1 mea- 
surements, and, based on our results, we suggest that the incidence of the magnetic field in various 
classes of stars should be revised. 
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1. INTRODUCTION 

The FORSl instrument (Appenzeller et al. 1998) of the ESO VLT has been extensively 
used for magnetic field measurements in various classes of stars. The large majority 
of FORSl measurements have been published in the literature, and we thought that 
a general catalogue would serve to obtain an overview (even though biased on target 
selection) of the incidence of the magnetic fields in various kinds of stars. However, a 
catalogue compiled using published material would suffer from the lack of homogeneity 
in the way data have been treated. Furthermore, over time, new ideas for data reduction 
and quality checks have improved the reliability of FORSl magnetic measurements, 
which calls for a revision of earlier data. We also noted that the literature of FORS 
magnetic field measurements includes a certain amount of controversial detections, e.g.: 
(i) FORS detections not confirmed by observations obtained with other instruments; (ii) 
FORS detections not confirmed by further FORS measurements; (iii) FORS detections 
not confirmed by the re-analysis of the same data by other groups; (iv) FORS detections 
not confirmed by the re-analysis of the same data performed by the same group who 



originally claimed detection. 

The release of the FORS pipeline for spectropolarimetric data (Izzo et al. 2010, 20 11) 
gave us the opportunity to develop a reduction method potentially more accurate than 
what available in the past. Compared to the reduction of raw data coming from individual 
observing runs, the mass-production of reduced spectra offers a few advantages in terms 
of quality of the final products, as the results can be checked on a large statistical basis. 
Therefore we have decided to tackle the task of re-reducing the entire FORS archive 
of spectropolarimetric data obtained in single target mode ("fast mode"), and compile a 
new catalogue. In this paper we report on some of our most important findings. 

In Sect. 2 and 3 we give details about data reduction and field measurements tech- 
niques, and in Sects. 4 and 5 we present our results and conclusions. 

2. MAGNETIC FIELD MEASUREMENTS 

The mean longitudinal magnetic field is obtained from the Stokes parameters I and V 
using the relationship 

V old/ 

y = -geffCzA2- — (5,), (1) 
where geff is the effective lande factor, and 

Cz=— (~4.67x 10-^^A-i G"^) (2) 

where e is the electron charge, me the electron mass, c the speed of light. It should 
be reminded that Eq. (1) is subject to several important limitations: (i) it is valid only 
in the limit of a magnetic field weak enough that Zeeman splitting is small compared 
to the local spectral line width (i.e., typically for fields of the order of a kG or less for 
optical metal spectral lines, or of order 10 kG for Balmer lines); (ii) it applies to isolated, 
unblended lines; (iii) the value of gefF varies by a significant amount from line to line; 
the use of an average value for all lines means that the actual (5^) value varies from the 
computed value in individual lines by up to ~ 25 %. 

Bagnulo et al. (2002) proposed that Eq. (1) could be applied to spectropolarimetric 
data obtained with FORSl, using a least-squares technique that was subsequently ex- 
tensively adopted in many surveys. Although field measurements obtained with Eq. (1) 
cannot be used for detailed modelling, this technique is certainly valuable to decide 
whether a star is magnetic or not, and to estimate the field strength with a reasonable 
approximation. 

3. FORSl ARCHIVE DATA REDUCTION 

All archive data obtained in fast mode have been pre-processed using the ESO FORS 
pipeline (Izzo et al. 2010, 201 1). Although this software tool is capable of delivering the 
final Stokes profiles, we preferred to develop a few algorithms that further improve the 
quality of the results. Therefore, instead of adopting the final products of the pipeline 
(i.e., the Stokes profiles), we have used the individual fluxes corresponding to the parallel 



and perpendicular beams split by the WoUaston prism (obtained at the various positions 
of the retarder waveplate), and combined them as explained below. 

We adopt the same formalism used in Bagnulo et al. (2009), i.e., /" and /-*- are 
the fluxes in the parallel and in the perpendicular beam of the polarization analyser, 
respectively, Py = V// is the circular polarization normalised to the intensity, and Ny 
is the null profile, a quantity that was introduced by Donati et al. (1997) , and that 
is representative of the noise of Py (see Bagnulo et al. 2009). As a quality check, a 
magnetic field can be deduced from the null profiles using the same formulas adopted to 
measure (B^. The diagnostic content of this null field value, {N^, will be discussed in 
Sect. 3.2. 

The extracted fluxes /" and /-*- are usually combined to obtain the Py and A'v profiles 
using the formulas of the difference methods given in Eqs. (A2) and (A7) of Bagnulo et 
al. (2009), which for convenience we reproduce below: 
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where aj belongs to the set {—45°, 135°}. We note that instead of setting the A/4 
retarder waveplate to all four possible angles, most of the observers preferred to set it 
only to the angles -45° and +45°. 

In practice, in many cases we found the Py profile slightly offset from zero, even 
when no circular polarization of the continuum was expected. A possible explanation 
is cross-talk from linear to circular polarization, as discussed by Bagnulo et al. (2009). 
However, we often found slight but noticeable offsets also in FORS data for stars that are 
not linearly polarized. These offsets can be naturally explained if the ratio between the 
transmission functions in the perpendicular beam, and the transmission function in the 
parallel beam, does not remain constant as the retarder waveplate is set to the different 
position angles. The ratio between the transmission functions can be estimated from the 
ratios of the /" and /-*- fluxes in each frame, and analytically propagated to the final 
expressions of the Py profiles using Eqs. (3) of Bagnulo et al. (201 1). 

Both Py and A^^ profiles show occasional spikes that occur in the same wavelength 
bin. Most of these spikes are probably produced by cosmic ray hits, and, if not removed, 
may lead to spurious detections of magnetic fields, or at least decrease the precision of its 
determination. As a remedy, we clip the Py profiles by discarding from the computation 
of the (B^) values those points for which the (rectified) Ny value departs from zero by 
more than 3 a. 

Figure 1 shows the results obtained for the FORS observations of the well known 
magnetic star HD 94660. 



3.1 Errors bars 



The uncertainties of the (5^) measurements can be calculated according to a well 
established error theory; assuming that each individual step of the measurement process 
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FIGURE 1. The polarized spectrum of the magnetic Ap star HD 94660. The top panel shows Stokes 
/ (black solid line, in arbitrary units, and not corrected for the instrument response), Py (red solid line 
centred about 0), and the null profile (blue solid line, offset by -0.75 % for display purpose). Light blue 
bars represent the error bars of Py. The slope of the red interpolating lines in the bottom panels gives the 
mean longitudinal field from Py (left panel) and from the null profile (right panel, the latter being expected 
zero), both calculated using the H Balmer lines only. (After Bagnulo et al. 201 1). 



can be linearised, one can propagate the error bars of the measured flux to the Py 
profile (see, e.g., Eqs. (A3-A6) of Bagnulo et al. 2009), and then obtain the a^g_^ error 
bar associated to the (B^) measurement. However, since there are indications that (5,) 
error bars may be underestimated, we have preferred to adopt the more conservative 
approach to multiply the analytical expression for C7^g„^ by the square root of the reduced 

X^- The validity of this approach has been tested with Monte Carlo simulations. We 
have considered all raw frames associated to a certain number of measurements, and 
calculated the {B^) and (A^,) values and their error bars. Then we have re-considered 
the same series of frames and scattered each pixel of each scientific frame according 
to a Gaussian distribution with a equal to the error pertaining to that pixel. We have 
re-reduced the so altered set of frames, and measured the {By) and (A^^) values from 
the H Balmer lines only and from the full spectrum. We have repeated this procedure 
5000 times, and obtained for each set of frames two distributions of (B^) values and two 




FIGURE 2. Monte Carlo simulations compared to predictions of the "standard" error theory for a 
strongly magnetic star (HD 94660) and a non magnetic star (HD 96441). Top panels refer to the (Bj) 
and (Nz) values obtained from the analysis of the H Balmer lines, and the bottom panels to the values 
obtained from the analysis of the full spectrum. Blue shaded histograms show the distributions of the {B-) 
and (Nz) values obtained after Monte Carlo simulations. Red solid lines show Gaussian curves centred on 
the original {B^) and (A^-) values with a calculated via usual standard error propagation. 



distributions of (A/,) values. We have compared these distributions with the Gaussian 
curves centred on the {B^) and (A^^) values measured from the unaltered frames, with o 
calculated via standard propagation theory. The results for two specific cases are shown 
in Fig. 2. The discrepancy observed for the case of (B^) obtained from the analysis of 
the full spectrum of HD 96441 (left bottom panel of Fig. 2) may be due to the fact that 
our error bar calculated analytically is weighted by a reduced X^- This value is definitely 
> 1 because Eq. (1) is not strictly valid for blended metal lines when the field is strong. 
All the remaining cases shown in Fig. 2 confirm that (B^) and (A^^) errors due to photon 
noise are correctly calculated. 




FIGURE 3. The actual distribution of the values (A^-) normalised to their error bars (shaded histogram) 
compared to their theoretical distribution, i.e., a Gaussian with (7=1. Vetical bars show the ±3 <J limits. 
(After Bagnulo et al. 201 1). 



3.2 Uncertainties not due to photon-noise 

There are indications that error bars due to photon-noise underestimate the actual 
uncertainties. 

As mentioned in Sect. 1, in the literature there are non negligible discrepancies be- 
tween the {B^ values obtained by different authors who analyse the same FORS dataset. 
Note also that sometimes, the same authors have revised their own field determinations. 
This is probably due to the fact that {B^ measurements depend, sometimes crucially, 
on the numerous specific choices taken during data reduction. For instance, the (5^) and 
(A^^) values depend on the algorithm adopted for spectral extraction, on how data are 
clipped to remove the effects of cosmic rays, on whether flat-field correction is applied 
or not. A complete description of all these issues would be outside of the scope of any- 
thing but a detailed monograph, hence here we limit ourselves to refer the reader to a 
few numerical examples given in Bagnulo et al. (201 1), and to warn that changes in the 
data reduction algorithm may lead to change of the {B^ values comparable to 1 a. 

A global inconsistency emerges from the analysis of the full FORS archive dataset, 
even when all data are reduced by consistently adopting the same algorithm, i.e., an 
unusually large number of field detections in the null profiles. Figure 3 shows the 
distribution of the null field values normalised to their error bars compared to the 
expected Gaussian with a = 1 . This anomaly could be explained by a number of reasons, 
e.g., (i) that for some dataset, data reduction needs an extremely refined tuning; (ii) stars 
change radial velocity during an observing series; (iii) during the observations there are 
instrument flexures. It is likely that all these issues play a role, but we should note that 
in particular, instrument flexures, if present, would be completely out of user control. 
A flexure as small as a quarter of a pixel cannot be detected with the monitoring of the 
position of sky lines, because most of the observations were obtained with very short 
exposure times, and because in general the observed spectral range includes only one 
telluric line. Yet, simple numerical simulations show that even such small flexures may 
change the (5^) and (A'^) values by an amount comparable to that of the error bars due 
to photon-noise. Occasional flexures of the order of 1 pixel could be responsible for the 



observed high frequency of (A^^) detections (about 3 %-5 %). 

4. RESULTS 

The magnetic field measurements obtained with FORSl have been directed at three main 
goals: detection of previously unknown fields, measurements of fields in stars already 
known to be magnetic, and establishing upper limits to fields in stars which are not 
magnetic. 

Our discussion of reduction procedures for FORSl spectropolarimetric data shows 
that using different reduction procedures to infer magnetic field strengths leads to differ- 
ent results. These include changes in both the value of the mean longitudinal field (B^) 
and the associated uncertainty a^g^^ obtained from any specific observation. 

This finding is probably not a serious concern for measurements of fields at the tens or 
hundreds of o level, as these measurements are much more sensitive to such factors as 
the lines available in the spectrum, the spectral region used, and the degree of chemical 
and inhomogeneity in the atmosphere. 

On the other hand, underestimating the actual measurement uncertainty can have quite 
important effects when the issue is whether or not a field has been detected at all in a 
particular star or class of stars. Failure to recognize that computed uncertainties may 
not fully represent the range of measurement values that may be taken by (5^) even 
in the absence of a real stellar field can lead to spurious discoveries. These in turn can 
lead to much effort being expended to analyse or explain the presence of fields that are 
not in fact there, or at least have not been reliably detected. It is in this context that 
the most harmful effects of incorrect measurements, or underestimated measurement 
uncertainties, can occur. 

Because of problems that spurious discovery announcements can cause in a research 
field, we have carried out a systematic analysis of new magnetic field detections obtained 
at less than about the 6o level. In the following we report on some of our findings. Our 
complete revision of FORSl magnetic data, which will include more classes of stars, 
will be published in two forthcoming papers (Bagnulo et al. 201 1, and in preparation). 

Herbig AeBe stars. The first measurement of a magnetic field in a Herbig Ae 
star obtained with the FORSl instrument was reported by Hubrig et al. (2004), who 
measured (B^) = —450 ±93 G in HD 139614. This discovery was deemed as spurious 
detection by Wade et al. (2005) (who from the same FORSl data measured (B^) = 
— 150±50G), but re-affirmed as real by Hubrig et al. (2006c). A new reduction by 
Hubrig et al. (2007a) led to the revised value of {B^) — — 112±36G, and our re- 
reduction of the discovery observation ((5^) = —84 ± 65 G) confirms that no field is 
detected in this measurement. The same star was observed several times using the high- 
resolution spectropolarimeter ESPaDOnS (Wade et al. 2005, 2007), and no magnetic 
field was ever detected. 

The results of a FORSl survey of 49 Herbig AeBe stars were described by Wade 
et al. (2005; 2007). Magnetic fields were reported in two stars: HD 72106A (whose 
companion HD 72106B is a Herbig star), and HD 101412. The uncertainties of the 
measurement of these two stars have increased in our re-reduction of these data, but both 



fields are still detected at 3 a level. These discoveries are corroborated by independent 
observations obtained with ESPaDOnS. Wade et al. (2005; 2007) report also marginal 
detections in other three Herbig stars, but these detections were never confirmed by 
subsequent observations, nor by our new data reduction. 

Further surveys of Herbig AeBe stars using FORSl were described by Hubrig et al. 
(2006c; 2007a), who reported the discovery of magnetic fields in two new Herbig stars. 
Hubrig et al. (2009b) report new magnetic field detections in six more Herbig stars 
(almost 30 % of the observed sample). In our re-reduction of all these data, only three 
measurements remain signficant. Detections in HD 144432 and in HD 144668 (at 3 o 
level) might be real, but require confirmation. A detection at the 5 — 6a level of a field 
in HD 150193 appears stronger in our re-reduction than originally published. However, 
a single ESPaDOnS observation with an uncertainty of about 15 G detects no field (E. 
Alecian, private communication). In conclusion, most of FORSl field discoveries of 
Herbig AeBe stars appear to be spurious, and detectable magnetic fields appear to occur 
in only < 10% of all Herbig AeBe stars observed with FORSl. 

HgMn stars. Some of the peculiar A stars of the HgMn type, and in particular a And 
= HD 358, have been repeatedly observed with FORSl. Based on these data, Hubrig et 
al. (2006b) have claimed field detections in four HgMn stars. 

Numerous investigations of magnetism in HgMn (e.g., Borra & Landstreet 1980, 
Shorlin et al. 2002, Makaganiuk et al. 2011) have found no convincing evidence for 
fields in this class of stars with uncertainties often in the range of a few G up to a few 
tens of G. Furthermore, a And has been specifically investigated with the MuSiCoS 
(13 observations, typical a ~ 30 — 60 G) and ESPaDOnS (5 observations, typical o ~ 
6 — 19 G) spectropolarimeters without any significant detection Wade et al. (2006). Our 
new reductions have increased the standard errors of all FORSl measurements by factors 
of 1.3, and none of our results shows a field at the 3 a level. 

j8 Cephei pulsators and Slowly Pulsating B stars. Hubrig et al. (2006a) and Hubrig 
et al. (2009a) have reported the discovery of a field in the j8 Cep star CMa = 
HD 46328, based on altogether 13 detections at approximately the 6 or la level. The 
field always appears to be close to +350 G. This detection is supported by a single field 
detection by Silvester et al. (2009) using ESPaDOnS at about the 30<7 level. Our re- 
reduction of these data also confirms the presence of a field. This discovery appears very 
robust, and clearly indicates the capability of FORSl to detect rather modest magnetic 
fields. 

Major FORS 1 surveys of both j8 Cep and SPB stars have been carried out by Hubrig et 
al. (2006a) and Hubrig et al. (2009a). Nearly 70 stars were observed for magnetic fields 
using the FORSl spectropolarimeter (19 known and suspected j8 Cep stars, 50 known 
and suspected SPBs). Fields were reported to have been detected in five j8 Cep stars 
and 26 SPB stars. If these results were confirmed they would be quite important. Such 
high incidence of detected fields would strongly suggest that magnetism is intrinsically 
connected with the pulsation phenomenon in early B stars, as it is in the cool rapidly 
oscillating Ap (roAp) stars. 

However, all reported field detections in j8 Cep stars, except those of CMa, have 
mostly become insignificant in the new reductions. Only two field measurements, one 



each for HD 74575 and HD 136504, are still barely significant at the 3(7 level. Of the 
more than 40 field detections reported in SPB stars by Hubrig et al. (2006a) and Hubrig 
et al. (2009a), all but five have decreased to non-detections in the new reductions. Our 
reduction suggest that fields might be present in HD 53921, HD 152511 (with two 
still apparently significant detections), and HD 208057. However, some or all of these 
detections could be manifestations of the occasional outliers found in the FORS 1 data. 
Note that a field has also been marginally detected by Silvester et al. (2009) in a single 
ESPaDOnS measurement of HD 208057. 

Based on the re-reductions of the FORSl data as well as other published material, 
we conclude that as in other kinds of upper main sequence stars, magnetic fields are 
relatively rare in /3 Cep and SPB stars. Furthermore, there is no significant case for 
considering that the pulsation properties of these B pulsators are intrinsically connected 
to the presence of weak magnetic fields. 

0-type stars. Hubrig et al. (2008) carried out the first large survey of O-type stars 
with FORS 1, and reported field measurements of 13 stars (including one star previously 
reported by Hubrig et al. 2007b). Field discoveries were claimed for five of these stars. 
Our reduction confirms a marginal detection in only one of these stars, the 06.5f?p star 
HD 148937. The field of this star has been detected in three more observations with 
FORS2 by Hubrig et al. (2011) and with ESPaDONs observations obtained by Wade 
et al. (2011). This discovery is clearly real, although the field of this star, always about 
—250 G, is close to the limit for reliable field detection by FORS 1 for such stars. 

Central stars of planetary nebulae. Jordan et al. (2005) carried out six field mea- 
surements of four central stars of planetary nebulae (CSPN), and concluded that at least 
two of these stars have highly significant fields in the kG range. This result supported 
the idea that such fields have an important shaping effect on the the planetary nebulae 
themselves, and suggest that there must have been important loss of magnetic flux be- 
tween this evolution state and the white dwarf state where stars with as much flux as 
that inferred for the CSPN are relatively rare. However, this result was called into ques- 
tion recently by Leone et al. (201 1), who obtained one new measurement of each of the 
CSPNs NGC 1360 and LSS 1362, detecting no significant field in either star. Leone et 
al. (2011) also re-reduced the older data of Jordan et al. (2005), finding no significant 
fields. Our new reductions are completely consistent with the re-reductions carried out 
by Leone et al. (2011). None of the original field measurements of Jordan et al. (2005) 
reveals a significant field. As remarked by Leone et al. (2011), there is now no signifi- 
cant evidence for coherent magnetic fields in the central stars of planetary nebulae, and 
the current best upper limits on {B^ are at roughly the 1 to 2 kG level. 

5. CONCLUSIONS 

FORS is an instrument perfectly capable of performing field measurements of the order 
of at least 250-300 G even in faint and rapidly rotating stars. However, for weaker fields, 
it is well possible that limits intrinsic to the technique, or even to the instrument itself, 
are reached. In particular, we have shown that (i) different data reduction techniques lead 



to results that differ for values comparable to (and in some cases larger than) photon- 
noise error bars; (ii) instrument flexures as small as a quarter of a pixel may affect the 
field measurements by an amount again comparable to a typical photon-noise error bar. 
Since FORS is mounted at the Cassegrain focus, such small flexures might occasionally 
occur. We conclude that, before being accepted as real, many recent claims of weak field 
detections obtained with FORS in various kinds of stars need to be confirmed by new 
and repeated measurements obtained both with FORS itself and, independently, with 
other instruments. 
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